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I. INTRODUCTION AND MOTIVATION 

Quantum computation (QC) can be implemented 
in various frameworks, such as the standard circuit 
model {Tl, adiabatic evolution [2], manipulation of exotic 
anyons in topological phases, and local measurement on 
certain entangled states [3]. In the measurement-based 
model of quantum computation (MBQC) j3H§], only cer¬ 
tain entangled states are known to provide the capability 
for driving a universal quantum computation via local 
measurement, such as the cluster state on the square 
lattice j7]. A complete classification of entanglement 
structure that enables MBQC remains a challenging open 
question. Moreover, whether these entangled states arise 
as unique ground states of short-range gapped Hamil¬ 
tonians is relevant for robust resource state creation and 
possibly further protection during computation [8] . From 
this latter viewpoint, cluster states, unfortunately, can¬ 
not be the unique ground state of two-body interacting 
qubit Hamiltonains B !, albeit they can be approxi¬ 
mately [TU] . The key obstacle for complete character¬ 
ization is that there is no simple physical observables 
(or order-parameter-like quantities) that is generic for 
answering whether a state is universal or not. Proving 
either universality or non-universality is thus in general 
highly nontrivial. In order to make progress toward com¬ 
plete characterization of universal resource states, it re¬ 
quires a substantial breakthrough in understanding the 
entanglement structure necessary for realizing universal 
gates. 

Until now the most complete characterization is for 
ID resource states, even though they are not universal 
for QC. This includes ID cluster state, ID spin-1 AKLT 


state, and matrix-product states of certain yet general 
forms mHHi- However, 2D and higher dimensions are 
much less understood. After the disovery of the clus¬ 
ter state on the square lattice, it was recognized that the 
generalization of the cluster state—the graph state—also 
provides universal resource on various other 2D regular 
lattices [IS] . Furthermore, it was also shown how to char¬ 
acterize quantum computational universality for graph 
states on faulty lattices [L6] as well random 2D planar 
graphs mi- The issue of the universality in the family 
of cluster or graph states is well understood due to their 
simple entanglement structure. Beyond this family of 
states, it seems only a handful of other entangled states 
are known to be universal mi (SHIS- Deciding whether 
a given quantum state is a universal resource for MBQC 
is still a challenge, let alone generalization to a family of 
states. No generally applicable strategies to answer this 
question are known. 

An interesting and emerging family that may po¬ 
tentially be as useful with regards to resourcefulness 
for MBQC is the so-called Affleck-Kennedy-Lieb-Tasaki 
(AKLT) states |2T|- Originally, the ID spin-1 AKLT 
chain provides the first evidence to support Haldane’s 
conjecture on the spectral gap property of integer-spin 
chains with rotational symmetry [22]. It also gives the 
first instance of the matrix product states [23] , and more 
recently serves as an example of a symmetry-protected 
topological ordered state [24] . The extension to two di¬ 
mensions, such as the ones on the honeycomb and square 
lattices, provides illustration of projected entangled pairs 
states [25] and models of spin rotational symmetry with 
potentially a finite gap above the ground state [21] . In 
three dimensions, some AKLT states can become Neel 
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ordered [2B]. Similar to graph states, AKLT states can 
be defined on any graph in terms of the valence-bond 
picture and local symmetrization, but the local spins can 
be of magnitude 1, 3/2, 2 or higher. Moreover, they are 
unique ground states of two-body interacting Hamilto¬ 
nians with suitable boundary conditions. Even though 
quantum computational universality in MBQC requires 
at least a 2D structure, the results on the ID spin-1 
AKLT state for theoretically and experimentally simu¬ 
lating one-qubit gates HU [271 E8] prompted the quest of 
universality in 2D AKLT states. However, the capabil¬ 
ity of full quantum computational universality was only 
established recently in the spin-3/2 AKLT state on the 
honeycome lattice [29] [30] and later on some other triva- 
lent lattices m , as well as some lattices that host spin-2 
and other lower spin (such as spin-3/2 or spin-1) hybrid 
AKLT states [32]. No AKLT states of uniform spin-2 en¬ 
tities have been known to provide universal resource, and 
on the contrary, the AKLT state on the kagorne lattice 
was argued to be non-universal. 

Here we demonstrate that the spin-2 AKLT state 
on the square lattice is indeed a universal resource for 
MBQC. This result adds a missing piece to a series of 
study mmm and gives rise to the following emerging 
picture that advances our understanding of the quantum 
computational universality in the valence-bond family. 
AKLT states involving spin-2 and other lower spin en¬ 
tities are universal if they reside on a two-dimensional 
frustration-free regular lattice with any combination of 
spin-2, spin-3/2, spin-1 and spin-1/2 (consistent with the 
lattice). Furthermore, geometric frustration can, but not 
necessarily, be a hinderance to the quantum computa¬ 
tional universality, and a frustrated lattice can be deco¬ 
rated (by adding additional spins) such that the resultant 
AKLT state is universal. Additionally, we argue that the 
spin-2 AKLT state on the three-dimensional diamond lat¬ 
tice is also a universal resource. The advantage of using 
a three-dimensonal resource state would be the possibil¬ 
ity of implementing fault-tolerant quantum computation 
with topological protection [55] , 

The family of AKLT states also provides the basis for 
generalization. They can be deformed so as to maintain 
universality in a range of the deformation parameter ^5] • 
We shall also argue that there is a finite region around 
the spin-2 AKLT point such that the ground state of a 
deformed AKLT Hamiltonian still supports a universal 
resource; see below in Sec. [V] Furthermore, frustrated 
AKLT states that are not likely universal (31] can be 
deformed in a such way that they are connected con¬ 
tinuously to a cluster state [35 . These can be used to 
study the connection of transitions in phases of matter 
and in quantum computational power da HH [351 [37]. 
With suitably chosen boundary conditions, AKLT states 
are unique ground states of certain two-body interacting 
Hamiltonians [21], some of which are believed to possess 
finite spectral gap [38, 39], including the spin-3/2 and 
the spin-2 AKLT states on the honeycomb and square 
lattices, respectively. A finite gap and the uniqueness of 


the ground state is a desirable feature for creating the 
resource state by cooling the system 0 [361 EQl HQ. 

The remainder of this paper is organized as follows. In 
Sec. [II] we describe the overall strategy for showing that 
an AKLT state is a universal resource for quantum com¬ 
putation. It consists of two steps, namely the mapping of 
the AKLT state to a random planar graph state by apply¬ 
ing a suitable POVM, and the numerical demonstration 
that, in the typical case, the resulting graph state can 
be mapped to a two-dimensional cluster state by further 
local measurements. These two steps are described in 
detail in Sections [ITT] and [IVJ respectively. In Sec. [V] we 
discuss one possible extension of our techniques to tran¬ 
sitions in quantum computational power away from the 
AKLT point. We conclude in Section [Vi] and also discuss 
possible future directions and potential experimental re¬ 
alizations. 


II. OVERALL STRATEGY 

Our goal is to show that any quantum computation 
that is efficiently implemented in the circuit model can 
also be implemented efficiently by a sequence of adaptive 
local measurements on a spin-2 AKLT state. In other 
words, we want to show that the spin-2 AKLT state is a 
universal resource for MBQC. 

The overall strategy for demonstrating this is the so- 
called quantum state reduction El: i.e., to show that, 
by local measurement, the state can be converted to a 
known resource state, such as cluster states, with finite 
probability, even in the limit of large system sizes. In 
our previous study of the spin-3/2 AKLT state [29], the 
reduction proceeds in three steps: 

1. Devise a pattern of local measurements that trans¬ 
forms the AKLT state into a graph state |G), where 
the graph G depends on the random measurement 
outcomes but is always planar. This proceeds in 
two sub-steps, namely (la) the creation of an en¬ 
coded graph state |G) by local generalized mea¬ 
surements, and (lb) a decoding of this graph state 
by local projective m easu rements. The support of 
each encoded qubit in |G) after step (la) is called a 
“domain” of the AKLT spin lattice C. The domains 
fluctuate in size. 

2. Show that a planar graph state |G) can be reduced 
to a 2D cluster state (the standard universal re¬ 
source), if the domains are all small and traversing 
paths through G exist. 

3. Numerically demonstrate that, for typical POVM 
outcomes in Step 1, the graph states produced sat¬ 
isfy the pre-conditions of Step 2. 

It is not a priori obvious that such procedure can work 
for the spin-2 or higher-spin AKLT states. Indeed it re¬ 
mains open whether one can even find a suitable POVM 
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for AKLT states with spin magnitude higher than two 
to reduce them to graph states. For the spin-2 case, the 
straightforward generalization of the POVM used in the 
case of spin 3/2 is no longer a POVM for spin 2. This 
is overcome by adding local POVM elements [32] such 
that the resulting measurement still maps to an encoded 
graph state, and the encoding can still be undone by local 
projective measurements. The additional local POVM 
elements amount to further projective measurement and 
hence disentanglement of the spins from the remaining. 
However, the new POVM does not guarantee that the 
resulting graph state |G) corresponds to a planar graph 
G, but the planarity of G is a requirement for Step 2 to 
work. Fortunately, the obstructions to planarity are lo¬ 
cal, and we can append a further round of measurements 
to remove them, thereby restoring planarity at the cost of 
reducing connectivity. We refer to this latter procedure 
as restoration of planarity by thinning ; see Sec. |III C| 

Once we obtain planar graph states, Step 2 of reduc¬ 
tion to 2D cluster states goes through unchanged for the 
spin-2 case. 2D Cluster states are the standard universal 
resource for MBQC, i.e., further local measurements on 
such a state can then implement any desired quantum 
circuit HHSj- 

Regarding Step 3, for correct numerical simulation 
of the measurement procedure, we require an efficient 
method for calculating the exact probability weights of 
the POVM outcomes. The probabilities for the outcomes 
F a and K a on each individual site are 4/15 and 1/15, re¬ 
spectively. Beyond those values, there are higher-order 
correlations between POVM outcomes on different sites 
which we cannot simply neglect. Even more seriously, 
some randomly assigned {F, K} do not occur, i.e., their 
probability is exactly zero. How do we know what out¬ 
comes can occur and with what probabilitity? It turns 
out that there is a closed-form expression for the exact 
probability weights which can be efficiently evaluated; 
see Sec. |IV A| for the expression and Appendix [C| for the 
proof. 

The most pronounced difference between the spin-3/2 
and spin-2 probability weights is that for spin 3/2 all 
possible combinations of POVM outcomes do indeed oc¬ 
cur with non-zero probability (except when the lattice is 
not bi-colorable, i.e., due to geometric frustration). This 
arises a consequence of the bi-colorability of the under¬ 
lying honeycomb lattice. For the spin-2 case, as already 
mentioned, certain combinations of POVM outcomes do 
not occur, i.e., have probability zero. The underlying 
spin lattice C (a square lattice) is still bi-colorable but 
this is no longer the deciding factor. 

Let us summarize the procedure to establish the uni¬ 
versality of the spin-2 AKLT state on the square lattice. 
The generalization of it to the diamond lattice is straight - 
foward but will not be carried out here. 

Procedure for establishing universality: 

1. Use the weight formula in Lemma 1 to sample typ¬ 
ical POVM outcomes. 
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FIG. 1. (a) AKLT state on the square lattice. Spin singlets 

| <j)) e = (|01) — 110))/v/2 of two virtual spins 1/2 are located 
on the edges of the square lattice, indicated by dashed lines. 
A projection at each lattice site onto the symmetric subspace 
of four virtual spins creates the AKLT state, (b) Illustration 
of domains and virtual qubits inside domains. Five of the 
domains are labeled: Cc,C m i,C M 2,C M 3,C/j4. 


2. Apply the thinning proceudre in Sec. |III C| to ob¬ 
tain assoicated random planar graphs. (One should 
not confuse the thinning procedure with the dele¬ 
tion used in checking the robustness of the graph 
connectivity, presented in the next step.) 

3. Check whether there is a traversing path and record 
the probability p span that this occurs. If so, exam¬ 
ine how robust the connectivity in the graphs. To 
do this, we employ the idea from percolation and 
delete every vertex with a probability Pdeiete an d 
record the probability p span that a traversing path 
still exists. 

If we can demonstrate, from the behavior of p span vs. 
Pdeiete for different sizes L , that there is a phase tran¬ 
sition (say at p* Ae i ete ), then according to the theory of 
percolation, the graphs that reside in the phase with 
Tdeiete < Tdeiete* (a.k.a. the supercritical phase) contain 
macroscopic number of traversing paths. As we shall 
demonstrate in Sec. ra that this is indeed the case, 
and hence the random graph states are universal, imply¬ 
ing the original AKLT state is also universal. 


III. REDUCTION FROM AKLT STATES TO 
GRAPH STATES 


Let us define the AKLT state on the square lattice. It is 
useful to view the spin-2 particle on each site is consisting 
of four virtual qubits. Each virtual qubit forms a singlet 
state, | (f) e = (|01) — 110))/v/2, with its corresponding 
virtual qubit on the neighboring site, with the singlets 
indicated by the dotted edges; see Fig. [l] In order to 
yield a local five-level spin-2 particle, a local projection 
on each site is made from the Hilbert space of four virtual 
qubits to their symmetric subspace, which is isomorphic 
to the spin-2 Hilbert space [21] ■ 
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A. Reduction from spin-2 entities to qubits: the 
generalized measurement 


The POVM we shall employ consists of three rank-two 
elements and three additional rank-one elements [52] : 

Fa = \/I (|5 “ = 2)(5q = 21 + l S a = -2)<S« = -2|)(la) 

Ka = ( lb ) 

where a = x,y,z and ]</>„) = (\S a = 2) ± \S a = —2))/a/ 2- 
The F”s are straightforward generalization from the spin- 
3/2 case [2S], but they do not give rise to the complete¬ 
ness relation, which is required for conservation of prob¬ 
abilities. By adding K’ s, it can be verified that the com¬ 
pleteness relation is satisfied: J2 a = 

1, i.e., proving that there only six possible outcomes as¬ 
sociated to F a - Xt y^ z and iv a =x, y,z- A spin-2 particle in 
state | <f>) that undergoes such a generalized measurement 
becomes either i 7 ’a|4 > ) or K a |$). 

The reduced density matrix for a single site of the 
AKLT state is a completely mixed state, and there¬ 
fore, each unwanted type occurs on average with a prob¬ 
ability 1/15. We note that the outcomes of K re¬ 
sult in projection to a one-dimensional subspace (in¬ 
stead of a two-dimensional subspace that could form 
the basis of a qubit), they are thus regarded as “un¬ 
desired” or “unwanted” outcomes. An unwanted out¬ 
come associated with K thus occurs with probability 
p err = 3 x 1/15 = 1/5. However, as we shall see be¬ 
low in Sec. lIV Al that not all POVM outcomes associated 
with sets of {F a r v \,Kp/ w )} occur with non-zero proba¬ 
bility, due to the correlation present in the AKLT state. 
Below we discuss the effect of F and K outcomes. 

For bi-colorable graphs, such as the square lattice, any 
all-F POVM outcome can occur. We have previously 
shown that in this case, the post-measurement state 

l G o({A})) ~ (§Z) F a ( v )\il> AKLT } (2) 

v€ C 

is an encoded graph state HUES], whose properties are 
described below. 


1. Each domain on C supports a single encoded qubit, 
i.e., the domains D C C are the sites or vertices of 
the graph Go, with the encoding as described in 
Tabl e [TJ The encoded qubits form a graph state 
|Go). When there is no confusion, we shall not 
distinguish between the graph state |Go) and its 
encoded version |Gq) and omit the labeling {F}. 


2. The graph G'o has an edge between the vertices 
v(D ) and v(D'), if the domains D and D' are con¬ 
nected by an odd number of edges in C. 


3. Be D a domain of type T £ {a", y, z} with n a neigh¬ 
bouring domains of type a. The stabilizer opera¬ 
tors for such a graph state are shown in Eq. (B1) in 


terms of encoded logical operators. They are char¬ 
acterized by the so-called stabilizer matrix, and in 
the case of graph state, is given via the adjacency 
matrix Ag 0 of the graph Gq. It is seen that when 


n y mod 2 = 1, for T = x, 

n x mod 2 = 1, for T = y, 

n y mod 2 = 1, for T = z, 

the stabilizer operator /C n has a logical Y operator 
at the support of D. This means that the graph 
Go has a self-loop attached to the domain D , i.e., 
(Ag 0 )d,d = 1- 

We recall the definition of a “domain”. A domain is 
a maximal set of neighbouring sites in the lattice C for 
which the outcome of the POVM Eq. (|T|) is F a or K a 
with the same a. That is, there are domains of x, y and 
x-type, and neighbouring domains must be of different 
type. The self-loop is a convenient picture to visualize 
the graph. But we can perform local logical rotation to 
transform Y to X so as to remove the self-loop, then 
the resulting stabilizer operators will be in the canon¬ 
ical form. (Such rotation will also change the basis of 
logical measurement.) Moreover, we shall often not dis¬ 
tinguish between an encoded X or Y operator from the 
corresponding X or Y operator, unless necessary. We 
also note that the stabilizer operator for graph states 
is usually defined as K-d = FXjj <S>D'enb(D) Zd 1 , where 
nb(D) denotes the set of D ’s neighbors, i.e., those ver¬ 
tices connected to D by an edge. We shall refer to such 
stabilizer operators as being in the canonical form and 
we have allowed additional signs in the definition. The 
graph state |G) can thus be defined by K,d\G) = |G) 
for all vertices D. However, under a local basis change, 
the stabilizer operator can be transformed to the form 
K-d = ±Fd <S>_D'enb(D) Zd 1 , where the operator at D is 
the logical Y. The above point 3 is to determine exactly 
what form of the stabilizer is for each domain belonging 
to the graph state |G 0 ); see also Eq. (Bl). 

We discuss the effect of K ’s in the next subsection. 
But let us remark that the effect of the POVM ([ij is to 
produce an encoded graph state |G) corresponding to a 
graph G with adjacency matrix Aq- In the previous case 
of spin 3/2 it was sufficient to identify the graph state 
only up to local unitary equivalence. For the spin-2 case, 
due to the more complicated POVM and weight formula, 
this is no longer the case. In particular, we need to keep 
track of the self-loops in G and the eigenvalues of the 
stabilizer generators for |G). It is useful to use the graph 
state | Gq) as a reference point for subsequent discussions. 


B. POVM outcomes K a : domain shrinking and 
logical Pauli measurements 

We now discuss the effect of K' s, which can be rewrit¬ 
ten as follows, 

K a = ^l/2\<t>-)(<t>-\F a = y/2j3K a F a . 


(3) 
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We can thus think of the POVM Eq. |l]) as a two-stage 
process: first the outcomes on all sites are F’s, and then 
a number of sites are flipped to K or equivalently a pro¬ 
jective measurement is done in the basis | <jy~) and the 
result \(j)~) is post-selected. 

We shall denote by [F. K] the POVM outcomes on 
all sites, by Jp C C the set of sites where the POVM 
outcome is of A-type, and by Jk = £\Jf the set of sites 
where POVM outcome is of A'-type. Upon obtaining 
{F,K} we can deduce the state |G) that the original 
AKLT state is transformed to, 

|G({X, A'})) = (££) I\ a ( u ) (^) -F a (v) IV’aklt) 

u£.Jk 

\Jk\ 

0 IC(u))(Co)l0 F a(«)IV’AKLT), 

uGJx vEC 

where the state is not normalized and the probability of 
the set of POVM outcomes {F,K} occurs is 

P({F,K}) = {G({F,K})\G({F,K})). (4) 

It was shown previously [22] that 


POVM outcome 

Z 

X 

y 

stabilizer generator 

\ . \ [i] U] 

AiAjUz F) z 5 

AiAjOx Ox 

XiXjOy^CTy^ 

X 

<8>i=l <Tx ] 

02LV 1 

0JLW 1 

z 

\ ._[*] 
Ax<T z 

\ ,„w 

Ai(Jx 

\ ,_H 

A% (Jy 


TABLE I. The dependence of stabilizers and encodings on 
the local POVM outcome. \C\ denotes the total number of 
sites contained in a domain C and iXj = 1..4|C| label any 
two distinct virtual qubits in the same domain C (as there are 
four vitural qubits in a site; see also Fig. E>)- The square 
lattice C is bi-partite and all sites can be divided into either 
A or B sublattice, V(C) = A U B, and A, = 1 if the virtual 
qubit i £ v £ A and Xi = — 1 if i £ v' £ B. This is due 
to the negative sign in the stabilizer generator for a singlet 
1(—crh ] cr]f ] )|= |0)yjnr a n edg e (i, j). The logical Y 
operator can be defined as Y = —iZX. 


the effect of the POVM element K z on |Go) is different: 
the encoded qubit living on that domain is projected into 
an eigenstate of X (with eigenvalue —1). In general the 
effect of all |^)(^| (associated with K a ) in multi-site 
domains C, including the correct normalization factor, is 
equivalent to a logical measurement of X operator, 



/ 1 \ |£| —|vj 

(^) F a{u ' j \ijj AKLT ) = c 0 |Go>, (5) 


where cq is an outcome-independent overall normaliza¬ 
tion, V is the set of domains, £ is the set of inter-domain 
edges (before the modulo-2 operation) and |Go) is prop¬ 
erly normalized to have unit norm |20| . For the encoding 
using virtual-qubit picture, see Table I] 

Summarizing the above discussion, we have 


|G({F,A})) = c 0 



\S\-\V\+\Jk\ 



|G 0 ({A})), 


( 6 ) 


where |Go({A})) is assumed to be properly normalized. 
Without the additional operators <g) u6Jif I 

the analysis of the computational universality would be 
the same as in the spin-3/2 case. It is these operators 
that complicate the situtation. However, as we shall see 
below their effect is not serious. 

First, as an example, consider a 2 -domain with two 
sites, one of which is measured in F z and the other in 
AT. As Table |T] shows, the effect of the measurement in 
AT is a mere shrinking of the domain from two sites to 
one . The graph G underlying the encoded graph state 
|G) remains the same, only the encoding changes on the 
domain in question. Strictly speaking, the normalization 
is reduced by a factor of y/2, due to the post-selection of 
only the outcome. 

As a second example, consider a 2 -domain comprising 
a single site and this site is affected by a K z . In this case, 


P c = [l + O c ]/2l^l = [1 + (—1)1^1 X c ]/2l^l, (7) 

where \V C \ denotes the number of sites in the domain 
and we have defined O c = (— l)' Vc ' X c ; see below in Ap¬ 
pendix [B] (In that section we also derive the form of the 
stabilizer for any domain, and it is seen that it is not 
necessary in the canonical graph-state form with A at a 
given vertex and Z’s at neighboring sites.) In the canoical 
graph-state basis (CGSB) this measurement corresponds 
to either the logical X or Y basis. If there is no self-loop 
on this domain, then the measurement in the CGSB is 
in the logical X basis and we shall refer to this domain 
as an X-measured domain. If there is a self-loop on this 
domain, we shall refer to this domain as a Y-measured 
domain as the measurement in the CGSB is in the logical 
Y basis, for which the graphical rule is to perform local 
complementation before removing the vertex (see below 
in Sec. |b|). 

In fact the above two examples give the complete ac¬ 
count of the effects caused by the POVM outcomes I\ a . 
We need to discuss each domain separately, and distin¬ 
guish two cases: (a) Fewer then all sites in an cc-domain 
are affected by the POVM outcome K a . Then, the do¬ 
main is simply shrunk, and the graph G is unaffected, 
(b) All sites in an a-domain are affected by the POVM 
outcome K a . Then, the encoded qubit residing on that 
domain is measured in the X-basis; see Eq. |7|. 

If the latter happens, in terms of CGSB, measurement 
can be either a logical X' or a logical Y measurement, 
and the state resulting from such measurement is again 
a graph state, and the new graph can be deduced from 
simple graph rules [123; see Fig. [3] for illustration for Y- 
measurement. 
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FIG. 2. (color online) Part of a random graph for domains 
(solid circles), (a) The square indicates an X-measured do¬ 
main and the hexagon indicates a Y-measured domain. In 
this example, the two measured domains are neighbors, and 
the effect on the graph will induce non-planarity. A simple 
approach is to apply active Z measurement on those domains 
(indicated by the diamonds) that enclose these connected X 
or Y-measured domains, similar to the game of go. (b) The 
upshot of the active Z measurements will remove these X/Y- 
measured domains as well as active Z-measured domains but 
will restore planarity. 


C. Restoring planarity by thinning 

From the discussions above, we understand that the 
AKLT state after POVM measurement on all sites is 
transformed to a graph state. However, the associated 
graph is generally not planar. We previously established 
simple criteria for computational universality of random 
planar graph states mm, namely their corresponding 
graphs need to have a traversing path, and the domains 
need to be microscopic. The latter requirement for do¬ 
mains to be microscopic was checked numerically in sev¬ 
eral trivalent lattices HZIESIEI] but can be argued to 
hold using percolation; see Section [TVB| 

The second step in the computational procedure, after 
the POVM Eq. ([!]), is therefore a further round of active 
measurements with the purpose of restoring planarity of 
the encoded graph state. The non-planarity was caused 
by the ](/)&)(fta \ outcomes on all sites in the domains of 
the graph state |Go). The thinning procedure degrades 
the graph state as a potentially universal resource, but 
it simplifies the universality proof. In other words, these 
measurements are performed for the sole purpose of sim¬ 
plifying the reasoning. 

Our strategy for restoring planarity is to first remove 
connected POVM “measured” (regardless of whether it is 
X- or Y-measured) domains by actively measuring their 
enclosing/neighboring domains in the logical Z basis, so 
as to remove these connected “measured” domains [42j . 
Pauli Y-measurements on graph states have the effect of 
removing the qubit in question from the graph state. On 
the level of the corresponding graph G, the given vertex 
along with all edges ending in it are removed |Uj. Hence, 
Y-measurements can be used to excise regions of a graph 
state. (We remark that the encoded ^-measurements 
can be implemented locally on the level of sites of £, 
as required. This is guaranteed by the coding arising in 
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FIG. 3. (color online) Graph transformation rules on Y 
measurements on a spin-2 domains, (a) illustrates the case 
where sites 3 and 6 belong to the same domain (indicated 
by the rectangle) and each of them has a K a outcome such 
that the effect on the domain is an effective logical Y mea¬ 
surement; others sites are assumed to be of different domain 
types (i.e. Fp^ a ). (b) gives an example where four distinct 
domains are connected to center spin-2 site and that there 
is an additional edge between domains 4 and 5. To further 
make the graph remain planar, logical Pauli Z measurements 
are performed, (c) shows the case where three distinct do¬ 
mains are connected to the spin-2 site. (This case can arise, 
e.g., as one of the neighboring domains was deleted in the sec¬ 
ond step of (a) or (b) associated with other spin-2 site.) (d) 
and (e) exemplify the cases of, respectively, two and one do¬ 
main connected to a spin-2 site. Note that the above list does 
not exhaust all possibilities (due to other possible connections 
between neighbors, albeit planar) but just serves to illustrate 
that Y measurement can be treated to maintain planarity. 


the given setting, c.f. Table |IJ) Thereby, we remove all 
X-measured domains and isolated multi-site (i.e. those 
with more than 2 sites) Y -measured domains by the same 
procedure. (We note that the X and Y are referred to 
in the canonical graph-state basis.) The non-planarity 
caused by these POVM “measured” domains is recovered 
quasi-locally; see Fig. [2] 

Then we proceed to deal with the remaining isolated 
Y-measured domains which contain either one single or 
two sites (which can have at most 6 neighboring sites 
and hence domains). The effect of Y-measured domains 
on the graph is to apply local complementation before 
removing the vertices corresponding to the Y-measured 
domains. If the Y-measured domain has three or fewer 
neighboring domains, the local complentation still pre¬ 
serves planarity. But when the number of neighbors is 
four or more, we then actively apply Z measurement on 
some of the neighboring domains (see Fig. [3]) to maintain 
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FIG. 4. (color online) Illustration of (1) [top left] { F , I\} configuration, (2) [top right] domain formation and graph construction 
(before treating (A'j’s), (3) [bottom left] treating domains which are A'-measured domains, and (4) [bottom right] treating 
Y-measured domains. Note that the measurement observables are referred to in the CGSB. The result is a plannar graph with 
fewer domains. In (2) domains are formed by connected sites with same type a = x,y,z (also color-coded) of either F a or K a . 
These sites are linked by thicker lines of same color. The edges corresponding to the graph for the graph state of domains. 
Square boxes indicate the effect of K’s on the domain is a logical X measurement (i.e. A-measured domain) whereas the 
hexagon boxes indicate the effect of K’s is a logical Y measurement (in the CGSB) or, in other words, Y-measured domains. 
For those domains that contain at least one F, there is no change of the graph due to I\ ’s. Note that for spins on the boundary 
of the lattice, we imagine, for convenience, that they are attached to spin-1/2 particles (not shown) so as to have four neighbors. 


local planarity of the graph. (In principle, any multiple- 
site Y-type flipped domains can be dealt with. But they 
appear with a lower probability.) In Fig. [4j we give an 
example from our simulation and show explictly the steps 
in obtaining a planar graph. In Fig. [5j we give the frac¬ 
tion of both the X and Y measured domains (due to the 
K POVMs) in the graph Go and the fraction of the addi¬ 
tionally Z measured domains needed to restore planarity. 

In the end we are left with a planar graph state, whose 
graph may or may not be percolated. If for large enough 
system and with finite nonzero probability, the graphs 
obtained after the above procedure are in the supercrit¬ 
ical phase, then the resultant graph states can be used 
for universal MBQC, implying the original AKLT state 
is universal as well. However, if the graphs are not in 
the supercritical phase, then it is inconclusive. Our sim¬ 


ulations indicate that we need to use L of order 80 or 
larger in order to show that the graphs are in the su¬ 
percritical phase with high probability such as 90%; see 
Fig. [6} Of course, more sophisticated procedure to deal 
with multiple-site Y-type domains as well as X-type do¬ 
mains can reduce the size L needed in the simulations, 
as more vertices will be preserved. But fortunately, the 
simple procedure described above turns out to be suffi¬ 
cient. 


To carry out the simulations, we still need to sample 
the configuration {F, K} accordin g to t he exact distri¬ 
bution p({F, K}) [29] , In Section IV A we describe the 
formula for it and the proof is relegated to Appendix [C| 
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FIG. 5. (color online) The fractions of measured domains vs. 
size L. The quantity ro (with “square” symbols) is the ratio of 
the number of both X and Y measured domains to the number 
of total domains in |Go), and rq (with “circle” symbols) is the 
ratio of the number of additionally Z measured domains (due 
to the thinning procedure) to the number of total domains. 
“Triangle” symbols represent the sum ro + rq. The fractions 
are approximately independent of the linear size L. 


FIG. 7. (color online) The probability of a traversing path 
Pspan vs. the probability of deleting each vertex (in the ran¬ 
dom planar graph after the thinning procedure) Pdeiete for 
L = 120 (“circle”), 140 (“square”), 160 (“cross”), 180 (“tri¬ 
angle”), where N = L 2 is the total number of sites. The lines 
are used to guide the eyes. The threshold of Pdeiete is ap¬ 
proximately 0.142(3). The crossing for these curves indicates 
that there is a percolation transition from the supercritical to 
subcritical phase in the thermodynamic limit. 
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valid measurement outcomes (as see below by the incom¬ 
patibility condition). This latter complication sets the 
spin-2 case apart from the spin-3/2 case (in addition to 
the POVM itself). Employing the exact sampling also 
enables us to estimate the probability (at least the lower 
bound) of obtaining a universal resource state from per¬ 
forming the reduction procedure. We note that as long 
as the reduction procedure gives a finite, nonzero suc¬ 
cess probability in the large system limit then the orig¬ 
inal state is still regarded as a universal resource state 
(though of probabilistic nature). The weight formula 
that we discuss below will enable the exact sampling and 
hence simulations using it will give final words on whether 
the spin-2 AKLT state on square lattice is a universal re¬ 
source for MBQC or not. 


FIG. 6. (color online) The probability of a traversing path, 
Pspan vs. the linear size L (with N = L 2 the total number of 
sites). As L increases p apan also increases. This is obtained 
with exact sampling. 


IV. EXACT WEIGHT FORMULA AND 
SIMULATION RESULTS 

The exact sampling is needed, as random assignment of 
F and K POVM outcomes does not correctly reflect the 
correlation that these outcomes must obey due to multi¬ 
partite entanglement in the AKLT state [33]. Moreover, 
many of randomly chosen assignment of F and K are not 


A. The weight formula 

Let us recapitulate the notations introduced in 
Sec. |III B| Consider a spin-2 AKLT state on a bi-colorable 
lattice C (generalization to non-bicolorable lattices is pos¬ 
sible), and POVM elements F a and K a (a = x,y,z). 
Denote by Jp C C the set of sites where the POVM out¬ 
come is of E-type and by Jk = £\Jf the set of sites 
where POVM outcome is of A'-type. Here additionally 
we denote by Dk the set of domains where the number 
of A-type POVM elements is equal to the total num¬ 
ber of sites in the domain. Denote {F,K} the set of 
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POVM outcomes corresponding to F^\ and and 

the probability for such occurrence is p({F, AT}). 

We have also introduced the graph state |Go) in 
Eq. and we will denote its stablizer group as <S(|Go)). 
The stablizer operators /C M ’s for |Go) with respect to the 
encoding in Table 0 can be either ±X C (g) MeNb ( c ) or 
±Y C ^gNbfc) ( see Appendix[B|, where Nb(c) denotes 
the set of neighbors of vertex c. As explained in Sec. |IIIB| 
the effect of AT-type POVM elements on a strict subset 
of sites in a domain only shrinks the size of a domain, 
whereas A"-type POVM measurement on all sites in a 
domain in Dk amounts to the measurement (on |Go)) 
of an encoded logical X. In the CGSB, all stabilizer 
operators are of the form X c (S^gNbfc) but then the 
effect of A'-type POVM elements in a domain inside Dk 
amounts to the measurement of an encoded observable 
either X or Y , depending on the existence of a self-loop. 
Let us label the set of all domains (i.e. vertices of the 
graph Go) by V, the set of all inter-domain edges in C 
by £ and the set of all edges of G 0 by E. Note that E is 
obtained from £ by a modulo-2 operation [25], Now we 
introduce a \V\ x \Dk\ binary-valued matrix H with its 
entries defined as follows, 

H^ = 0, if[AC M ,O„]=0, (8a) 

H ltv = l, if{/C M ,O„} = 0, (8b) 

where is the stabilizer operator associated with the 
vertex (or domain) /i G V of the graph G 0 , 0„ is the 
operator defined in Eq. ([7]), and v G Dk ; see also Ap¬ 
pendix [b 1 Let dim(ker(fL)j denote the dimension of the 
kernel of matrix H. The utility of H is in the following 
lemma. 

Lemma 1 If there exists a set Q C Dk such that 
— Op € <5(|Go)), then p({F, A'}) = 0. Otherwise, 

/-. \ |£|-|V+2|./x|—dim(ker(ff)) 

p({F,7F})=c(-) , (9) 

where c is a constant. 

We subsequently refer to the above condition for 
p({F, A'}) = 0 as the incompatability condition. The in¬ 
compatibility condition implies that not all POVM out¬ 
comes labeled by F* and I\ a can occur. When there is no 
K outcome, Eq. En reduces to p = c 2 l y H £ l of previous 
results [12). The correlation of F’s and K’s at different 
sites is reflected either in the incompatibility condition (if 
it is met) or else in the factor dim(ker(Lf)). The proba¬ 
bility distribution of {F, AT} is thus very far from being 
independent and random. For the proof of the lemma, 
see Appendix [C| 

With the weight formula we can sample the exact dis¬ 
tribution of physically allowed POVM outcomes {F, AT} 
and carry out the procedure to restore planarity of the 
random graphs associated with the post-POVM states. 
To show that these graph states are universal, we need 


to show that (i) the domains are not of macroscopic size 
and (ii) these graphs reside in the supercritical phase 
(as the system size increases). For (i), we remark that 
the largest domain size can only be logarithmic with the 
system size N. This was confirmed earlier in our simu¬ 
lations of the honeycomb case. In fact, it is well-known 
from percolation that below the percolation threshold, 
the largest connected structure is logarithmic in N (by 
a previous study of percolation by Bazant [S]). This 
can be applied to our present study. If we only sample 
F’s, there are three types of outcome x, y , z and locally 
with the probability 1/3. Domains consist of connected 
sites of same type. For each type, it is a site percolation 
problem, at the occupation probability 1/3, smaller than 
the site percolation threshold 0.59 (cf. for honeycomb, 
the threshold is 0.69). This means that each type of do¬ 
main can be at largest logarithmic in N. Because 1/3 
is very far from 0.59 and the AKLT state has a small 
and finite correlation length, inclusion of correlation will 
not change the scaling. Furthermore, the K's can only 
decrease the domain size, so the relevant largest domain 
size is for sampling F’s only. 

The sampling is obtained by using the standard 
Metropolis algorithm for updating {F, A'} configurations. 
One notable distinction is that we need to check whether 
the next configuration satisfies the incompatibility condi¬ 
tion. If it is satisfied, we then abondon that configuration 
and generate another one until the incompability condi¬ 
tion is not satisfied. Then the configuration is accepted 
using the standard probability ratio, as was done in e.g. 
the spin-3/2 AKLT case [22 1. 


B. Numerical simulations 

Here, we first take typical random planar graphs after 
the thinning procedure and check the probability p span 
that a traversing path exists. We find that p span in¬ 
creases as the linear size L of the square lattice increases 
and approaches to unity eventually; see Fig. [6] This sug¬ 
gests that for L large enough, the random graphs re¬ 
sulting from the thinning procedure are percolated. To 
confirm this and examine the connectivity of the typical 
random graphs and perform site-percolation simulations, 
by removing every vertex in the graphs with a probabil¬ 
ity ^delete- It is important to note that we are interested 
in the connectivity for the random planar graphs result¬ 
ing from the thinning procedure, i.e. those graphs at 
Adeiete = 0, but as a means to characterize the robust¬ 
ness in the connectivity we examine how much we have 
to delete the vertices in order to remove all the traversing 
paths that were there and then use the obtained thresh¬ 
old value pjeiete as th e quantification. We expect that 
when ^delete is sufficiently large, there will not be any 
traversing path left. Indeed this is what we observe. 
Additionally, the crossing of curves in Fig. [7] for differ¬ 
ent sizes indicates that there is a percolation transition 
(at Pdeiete ~ 0.142) from the supercritical to subcritical 
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C. AKLT state on the diamond lattice 



(Pdelete-Pdelete*)L 1/nU 


FIG. 8. (color online) Data collapse in p spa n vs. (pdeiete — 
Pdeiete) ' L 1t/,v for different L = 120 (“circle”), 140 (“square”), 
160 (“cross”), 180 (“triangle”), where N = L 2 the total num¬ 
ber of sites. The lines are used to guide the eyes. We note that 
Pspan represents the probability of a traversing path and Pdeiete 
the probability of deleting each vertex in the random planar 
graph after the thinning procedure. The Pdeiete ~ 0.142 is 
extracted from the results in Fig. [7] and v = 4/3 is used from 
2D percolation. There is no tuning parameter. Data collapse 
is clearly seen (except for small deviation in the data from 
L = 120, possibly due to the small size) and this confirms 
that there is a continuous phase transition. 


phase in the thermodynamic limit. This then justifies the 
separation of two phases: (a) supercritical phase of perco¬ 
lation when pdeiete < Pdeiete j where there are macroscopic 
number of traversing paths; and (b) subcritical phase of 
percolation when Pdeiete > Pdeiete’ where no traversing 
path exists. Graphs residing in the supercritical phase of 
percolation contain a macoscopic number of traversing 
paths. This shows that our system sitting at Pdeiete = 0 
can be used to generate a network of entanglement that is 
universal for measurement-based quantum computation. 

As a further confirmation of the transition, we can 
rescale the horizontal axis via (pdeiete — Pdeiete)^ 1 ^ an d 
collapse all data points approximately on a universal 
curve; see Fig. [8] The exponent v = 4/3 is used for the 
correlation-length critical exponent of the 2D percolation 
universality class. There is no tuning parameter used and 
the data collapse demonstrates that the phase transition 
is continuous. By establishing the transition and the crit¬ 
ical point, we have thus showed that the random graph 
states we generated possess graphs residing in the su¬ 
percritical phase of percolation and hence these random 
graph states are universal for MBQC El- Therefore, the 
spin-2 AKLT state, from which the random graph states 
are generated by local measurement, is itself a universal 
resource. 


The technique for spin-2 square lattice can be applied 
to other planar lattices with the coordination number 
2 = 4 and be extended to 3D ones, such as the dia¬ 
mond lattice. For the latter lattice one can generalize 
the consideration to three-dimensional percolation and 
consider conversion of a 3D random graph state (in the 
supercritical phase) to a 3D cluster state, which can then 
be useful for providing a fault-tolerant implementation 
of MBQC [55] , Even though we do not pursue techni¬ 
cal demonstration in the present paper, we believe that 
the AKLT state on the diamond lattice is universal for 
MBQC. 

First of all, as the site percolation threshold of the dia¬ 
mond lattice is approximately 0.430, larger than 1/3, the 
largest domain size will be at most logarithmic in the to¬ 
tal number of the spins, and hence is not macroscopic. 
Secondly, the planarity that was the key obstacle in 2D is 
now relaxed to three-dimensionality. Thus fewer domains 
need to be actively measured in the canonical Z basis and 
the connectivity can be more easily preserved. The re¬ 
sulting 3D random graph states are thus expected to pos¬ 
sess sufficient connectivity in terms of percolation. From 
there, the reduction of 3D random graphs to 3D regular 
graphs could be carried out by local measurement. The 
sufficient connectivity implies that measurement-based 
quantum computation can be realized. The main advan¬ 
tage of using such a 3D resource state is that the quantum 
computation can be implemented in a topologically pro¬ 
tected manner and can tolerate high error rates [33:. We 
remark that there exist 2 = 3 regular structures in 3D 
with higher percolation thresholds than 1 /3 [fS] . AKLT 
states on these lattices also support fault-tolerant quan¬ 
tum computation, but these 3D lattices with 2 = 3 are 
not as natural as the diamond lattice. 


V. DEFORMED AKLT AND TRANSITION IN 
QUANTUM COMPUTATIONAL POWER 

In this section we discuss how our techniques and re¬ 
sults can be extended to inquire the question whether 
there is an extend region in the Hilbert space surround¬ 
ing the AKLT point such that the states can also provide 
universal resource. Is there a transition in terms of quan¬ 
tum computational power? In the following we shall pro¬ 
vide argument to support the affimative answers without 
carrying out detailed numerical simulations. 

Let us define a deformation operator 

T>(a) = i(|5 z = 2)<5 z =2| + \S z = -2)(S z = -2\) + (10) 
\S X = 1)(S X = 1\ + \S X = -1)(S X = -1\ + |S* = 0)<S z = 0[, 

where we shall restrict to the case a > 1/ \/3. The opera¬ 
tor V(a) is then used to deform the original spin-2 AKLT 
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Hamiltonian into 

H D (a) = Y / [V(a) i ®V(a) j ]h?;f hT [V{a) i ®V(a) j ], (11) 
(iJ) 


where /ifj <LT ’s are the two-body terms between pairs of 
neighboring sites (i, j) in the original AKLT Hamiltonian. 
Similar deformation was discussed in the spin-1 and spin- 
3/2 cases pSl B51 B3] . I' is easy t° see that the following 
one-parameter deformed AKLT state 


I Mo)) = [V~\a)f N I^klt) (12) 


is the ground state of the deformed Hamiltonian 
and it reduces to the original AKLT state at a = 1. 
Moreover, when a 1, | ipu(a)) is essentially the Neel 
ordered state, or more precisely, a superposition of two 
different Neel patterns, \S Z = 3/2) ® \S Z = —3/2) (g) • • • 
and \S Z = —3/2) g |5 Z = 3/2) g • • •, which does not have 
sufficient entanglement for MBQC, similar to the spin- 
3/2 case jS5|. As the spin-2 AKLT state on the square 
lattice (and the diamond lattice) is not Neel ordered, the 
wavefunction in (10) must possess a transition point in 


terms of phases of matter at certain value of a. 

The questions that concerns us here is whether there 
is a finite region around the AKLT point a = 1 such that 
the wavefunction (10) is still universal for MBQC. (An¬ 
other interesting question is whether the disappearance 
of the universality coincides with the valence-bond-solid 
to Neel transition, but we leave this for future investiga¬ 
tion.) We shall argue below that indeed there is a finite 
region around the AKLT point such that the ground state 
is universal. 

First of all, due to the deformation, the POVM in ( T|) 
cannot be directly applied but can be modified to work, 


F x (a) 

= F x D(a), F v (a ) = F y V{a), 

(13a) 


1 3 — a 2 


F z{a) 

- V 2a2 

(13b) 

K a (a) 

= K a v(a), 

(13c) 


with which one can verify that ]T) F ^(a)F a (a) + 
J2 a K^(a)K a (a) = 1. The physical intuition of how the 
above POVM can work for the purpose of MBQC is that 
it can be regarded as a two-step process: (i) first, the 
deformation T>(a) undoes the action T>~ 1 (a ) in |^d(o)) 
and converts it back to the AKLT state; (ii) the original 
POVM {F ai Kp} then acts on the AKLT state as before. 

Furthermore, the weight formula in Lemma 1 is then 
modified to 


p{{F a ,Kp}) = c(a) 


3 — a z 
2a 2 


p({F,K}), (14) 


where c(a) is an a-dependent normalization that is inde¬ 
pendent of {F, K}, riF z denotes the total number of sites 
where the POVM outcome is F z . and p({F,K}) is given 
by Eq. ^ and the incompatibility condition remains the 
same. The only difference in the weight formula is the 


additional factor involving uf z and when a = 1 the for¬ 
mula reduces to Eq. §• As we have found a finite per¬ 
colation threshold p)j elete for a = 1 case and the weight 
formula is continuous in a, we thus expect to have in¬ 
dependent threshold delete (°)> exce pt when the largest 
domain size becomes macroscopic. Thus there must exist 
a finite range around a = 1 such that the pj elete is nonzero 
and the largest domain size is microscopic. Therefore, 
the deformed AKLT state (10) is universal for MBQC 
in this region, and there exists a transition of quantum 
computational power as a increases, from being universal 
to non-universal. This argument also applies to the three 
dimensional case, and in particular, the diamond lattice. 

Even though we do not carry out simulations here to 
pin down the exact transition, we describe what may 
be achieved with present techniques. By studying the 
dependence of the largest domain size on the parame¬ 
ter a , we can extract at least the upper bound on the 
transition, a up , when the largest domain size begins to 
become macroscopic. Via the dependence of J?d e i ete (a) 
and the location, ai ower , at which Pd e i ete ( a iower) becomes 
zero, we can extract the lower bound of the transi¬ 
tion. The exact transition a t ra ns of the quantum com¬ 
putational power will then lie between the two bounds: 
Slower < dtrans < u U p- One could also determine the spon¬ 
taneous staggered magnetization to extract the transi¬ 
tion of the phase of the matter. It would be interesting 
to see how the two transitions differ. 


VI. CONCLUDING REMARKS 

The family of Affleck-Kennedy-Lieb-Tasaki states pro¬ 
vides a versatile playground for universal quantum com¬ 
putation. The merit of these states is that by appro¬ 
priately choosing boundary conditions they are unique 
ground states of two-body interacting Hamiltonians, pos¬ 
sibly with a spectral gap above the ground states. Here 
we have overcome several obstacles and shown that the 
spin-2 AKLT state on the square lattice is also a uni¬ 
versal resource for measurement-based quantum compu¬ 
tation. In particular, we were able to derive an exact 
weight formula for any given POVM outcome. Com¬ 
bined with a thinning procedure to restore planarity of 
random graph states, we performed Monte Carlo simula¬ 
tions and demonstrated that the assoicated planar ran¬ 
dom graphs from the procedure possess sufficient con¬ 
nectivity and reside in the supercritical phase. In partic¬ 
ular, our numerical site-percolation simulations showed 
that as the deletion probability increases (i.e., as the oc¬ 
cupation probability decreases) the system of the above 
random graphs makes a continuous phase transition from 
the supercritical phase of percolation to the subcritical 
phase of percolation. Moreover, the continuous phase 
transition is consistent with the universality class of the 
2D percolation. These demonstrate that typical random 
graph states obtained via our local measurement proce¬ 
dure on the AKLT state are universal for measurement- 
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based quantum computation. Thus, the spin-2 AKLT 
state on its own is also universal. 

One of the important enabling points for our proof is 
the spin-2 POVM. This POVM was used previously by 
us in considering hybrid AKLT states with isolated spin- 
2 and lower-spin entities [32]. Here, we were able to deal 
with the case that spin-2 particles are neighbors. An¬ 
other enabling point is the exact weight formula that we 
derived for arbitrary measurement outcomes according 
to a spin-2 POVM on all spins. This formula can be 
extended to Pauli measurements on stablizer states [13] . 
which may be useful for classical simulations of certain 
gates. Our weight formula also demonstrates the most 
pronounced difference between the spin-3/2 and spin-2 
cases: for spin-3/2 AKLT on bi-partite lattices, all pos¬ 
sible combinations of POVM outcomes do indeed occur 
with non-zero probability, whereas, for the spin-2 case, 
certain combinations of POVM outcomes do not occur, 
i.e., have probability zero. This shows that the POVM 
outcome in spin-2 case is very different from being ran¬ 
dom and independent, and the simulations from the ex¬ 
act weight formula are crucial in answering whether the 
spin-2 AKLT state is universal for MBQC or not. 

The emerging picture from our series of study on 
the quantum computational universality in the two- 
dimensional AKLT valence-bond family is as follows. 
AKLT states involving spin-2 and other lower spin en¬ 
tities are universal if they reside on a two-dimensional 
frustration-free regular lattice with any combination of 
spin-2, spin-3/2, spin-1 and spin-1/2 (consistent with the 
lattice). Additionally, the effect of frustrated lattice may 
not be serious and can always be decorated (by adding 
additional spins) such that the resultant AKLT state is 
universal. We conjecture that the result hold in three 
dimensions as well. 

We remark that an alternative approach to demon¬ 
strate the quantum computational universality is to ex¬ 
plicitly construct universal gates and show that they can 
be realized to simulate arbitrary quantum circuits. This 
can be done via the so-called correlation-space approach 
of MBQC as demonstrated in the spin-3/2 case [30] . 
One could carry out this procedure for the spin-2 case. 
However, similar to the spin-3/2 case, the key enabling 
ingredient is the POVM |l]), as well as the thinning pro¬ 
cedure in Sec. |III C| After these steps, the state becomes 
a graph state, and thus the gate construction in the corre¬ 
lation space can be implemented accordingly. Moreover, 
the question of whether there are sufficient gates that 
can be used to simulate arbitrary quantum circuits may 
still rely on the percolation argument and the numerical 
simulations done here. However, it may be possible to by¬ 
pass the thinning procedure and directly design gates to 
ensure that arbitrary quantum circuit can be simulated. 
But here we do not pursue this alternative approach. 

Simple states that are ferromagnetic or Neel ordered 
are believed not to possess entanglement structure use¬ 
ful for MBQC. AKLT states in 2D were shown to be 
disordered, i.e., not Neel ordered, and some (not includ¬ 


ing the diamond lattice) in 3D were Neel ordered f2Ej . 
Could it be that being disordered (and non-frustrated) 
is a sufficient condition for these AKLT states being uni¬ 
versal? We believe it is not. In terms of spin magni¬ 
tude, one expects that the larger the spin magnitude the 
more classical the system becomes, i.e., the effect of non¬ 
commutativity of spin operators becomes less and less 
important. For larger spin magnitdue S, despite being 
disordered, the system is becoming more classical. On 
this ground we suspect that for large enough S , AKLT 
states will eventually cease to be universal for MBQC. 
But what is the boundary of being universal and non- 
universal and what would be the decisive physical prop¬ 
erties for determining such a boundary? To address this 
question requires further insight. Nevertheless, our re¬ 
sults in the present paper show that the S = 2 case is 
still in the universal regime. 

One direction of generalization is to investigate the ro¬ 
bustness of the resource under small perturbations, e.g., 
slightly away from the AKLT Hamiltonian. We have 
discussed a particular deformation of the spin-2 AKLT 
Hamiltonian and the transition in quantum computa¬ 
tional power in Sec. [V] One could also consider defor¬ 
mations that take the spin-2 kagome AKLT state to an 
effective cluster state, turning a non-universal state to a 
universal one. It would also be interesting to consider 
the effect of finite temperature and at how high the tem¬ 
perature the system considered here can still support a 
universal resource ESZj- This is beyond the scope of the 
present paper and is left for future investigation. An¬ 
other direction one can study is the connection of the 
resourcefulness of the AKLT states (and their general¬ 
ization) to the symmetry-protected topological (SPT) or¬ 
der, the connection of which was recently found in one 
dimension 49.] ■ AKLT states serve as concrete exam¬ 
ples of SPT ordered states. However, what symmetry is 
needed and in what SPT phases can there be naturally 
protected universal gates? Progress along this direction 
in two and higher dimensions can potentially advance 
our understanding towards complete characterization of 
universal resource states. We also leave this intriguing 
question for future investigation. 

We wish to conclude by commenting on potential ex¬ 
perimental implementations. Bosons with multiple hy- 
perfine states can be used for high spins, and thus they 
are a natural candidate to realize possible AKLT-like 
Hamiltonians when placed in optical lattices [SEElC]- The 
difficulties lie in (1) tuning interactions to the desired 
Hamiltonian regime and (2) local controllability of single 
atoms. Significant experimental progress has been made 
in (2) that it is now possible to detect and image single 
atoms [52j [53]. One advantage of utilizing cold atoms in 
optical lattices is the ability to scale up the system. In 
addition to such a top-down approach, it is also possi¬ 
ble to build the resource states from bottom up, such as 
with entangled photons. The first experimental demon¬ 
stration of implementing gates in the ID AKLT state 
was carried out in the photonic system |28| . The key 
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idea relies on the mapping of the general AKLT state to 
a bosonic state (which has a Laughlin-like wavefunction 
in the coherent-state representation) [53] , in that 

I^AKLT (C,M)) ~ - b\a]) M \0), (15) 

(i,j) 


where |0) is the vacuum state, a\ creates a boson of type 
a (such as the horizontal polarization H of a photon) 
at site i, b t creates a boson of type b at site j (such as 
the vertically V polarizaed photon), and M denotes the 
number of singlets, with M = 1 being the original AKLT 
state. Therefore, creating an AKLT state is equivalent 
to placing singlet pairs of photons, such as | HV) — | VH), 
according to the valence-bond construction [2T] and the 
symmetrization of the photons are automatic [551 . In the 
ID AKLT state, there are two photons “per site” (or per 
mode), and the measurement of the effective spin-1 de¬ 
gree of freedom can be carried out if the two photons are 
indistinguishable [53 HZ]. There are three and four pho¬ 
tons per site on the 2D honeycomb and square lattices, 
respectively. Measurement on multiple indistinguishable 
photons that is equivalent to measurement in the spin 
basis can be carried out similarly [53 [57]. It is thus 
necessary to make those photons on the same site in¬ 
distinguishable, i.e., to achieve good mode matching. A 
proposal was made in realzing the spin-3/2 AKLT state 
on the honeycomb lattice and in implementing the key 
POVM [55]- Its generalization to the spin-2 case on the 
squre and diamond lattices is possible. The advantage of 
using entangled photons is that small-scale implementa¬ 
tion is already within the reach of current technology, but 
the disadvantage would be to scale up to a large system. 
Beyond the atomic-molecular-optical schemes, it is very 
recently proposed to realize AKLT and general valence- 
bond states in solid-state systems with, e.g., t 2g electrons 
in Mott insulators Ell- 
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Appendix A: The POVM expressed in terms of 
four-qubit representation 


spin-2 POVM are 


Fs = vf (l ° 04)(o04|+|1 * 4)(1 * 4|) 

(Ala) 

Fx=]/l(\+ m )(+n+\- m )(-n) 

(Alb) 

^ = vl {l ^ 4)(i04|+IH)04)(H)04|) 

(Ale) 

K x = \j\ |GHZ;)(GHZJ| 

(Aid) 

K x = ^j\ |GHZ-)(GHZ-| 

(Ale) 

K y = J ||GHZ;)(GHZ-|, 

(Alf) 


where |i /j® 4 ) is a short-hand notation for | ip,ip,ip,ip), 
equivalent to an eigenstate |S = 2 ,S a ) of the spin-2 op¬ 
erator S a with an eigenvalue S a = ±2 in either a = x, y, 
or z direction. Note that strictly speaking we should use 
different notations for these POVM elements, but since 
they are equivalent to those in Eq. (fT|), we retain the same 
notations. Moreover, cr z |0/l} = ±|0/1), a x \±) = ±|±), 
and a v \ ± i) = ±| ± i). The first three elements in 
Eq. (Al) are similar to those in spin-3/2 sites, except 
the number of virtual qubits being four, and correspond 
to good outcomes of type x, y and z, respectively. Asso¬ 
ciated with the last three elements, |GHZJ) = (|0000) — 
|H11 ))/V2, |GHZ-) = (|+ + + +>-| —--»/\/2, and 
\GHZ y ) = (|i, i, i, i) — | — i, -4, -4, -4))/y/2 are the corre¬ 
sponding states and they will be regarded as unwanted 
outcomes of type x, y, and z, respectively. But these 
GHZ states are at least eigenstates for certain product 
combination of Pauli operators. It can be verified that 
Ea F l F a + Eq K l K a = F s, where P s is the projection 
onto the symmetric subspace of four qubits, i.e., identity 
in the spin-2 Hilbert space. Let us also note that the I \’s 
operators can be rewritten as 


K a = v / i72|GHZ-)(GHZ-|Eq. (A2) 


Furthermore, there is a useful relation: 


(! 


[u;l] [u;2] [u;3] [u;4K 


IK 


|GHZ-)(GHZ"| =n ( 

(A3) 

where n a is a projection to a two-dimensional sub¬ 


space and is an identity operator on the code subspace 
(for the corresponding POVM outcome); specificially, 

n* = | + + + +)(+ + + + | + |-)(-|, 

n y = \i,i,i,i)(i,i,i,i\ + I - i,-i,-i,-i)(-i,-i,-i,-i\ 
and U z = |0000)(0000| + 11111)<11111 . The label b a de¬ 
notes the corresponding type b if a c = a ; see Table [jT] 


Appendix B: The exact form of stabilizer generators 


Expressed in terms of the four virtual qubits repre¬ 
senting a spin-2 particle, the elements that comprise the 


In this section we give the explict form of the stabilizer 
operator K-c for the domain labeled by C. It includes all 
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subtle plus and minus signs. The result is general for all 
states |Go) ~ |^aklt)j where F’s can be of 

arbitrary spins. This was already considered in the case 
of the spin-3/2 AKLT state [25], but the argument used 
there applies more generally. 

Let us first explain the notation. Consider a central 
vertex G € V(Go({E})) and all its neighboring vertices 
C/j, £ V'(Go). Denote the POVM outcome for all £-sites 
v £ G, G M by a c and a M , respectively. Denote by E y the 
set of /hedges that run between G and C^. Denote by 
E c the set of G-edges internal to C. Denote by V c the 
set of all qubits in G, and by V y the set of all qubits in 
C/j,. (Recall that there are 4 qubit locations per G-vertex 
v £ G, G m .) Extending Eq. (33) of Ref. |T7] to the spin-2 
case, we have 


ICc = 


(u(e))(v@) 


(-1 )ar K ~"cr. 


(- 1 ), 


(t>i(eO) (^(eO) 


A 4 eGEp, 

= i e p\ i 


e'£E c 


rr(n(e)) (v(e)) 
a u u a u 


A 4 

(«i(eO) _(« 2 (e 0 ) 
a b a b 


e' €E c 


We shall take the following convention for b as shown 
in Table [Tl] For POVM outcome a c = z , we take b = x; 
for a c = x, we take b = z; for a c = y, we take b = z. 
With this choice we have 


K, c = (-1)1®-!+^ 1^1 0(® eeEM A u(e) )Zlf-l 


a v (e) _v( e ) x 


e£E u 


It is convenient to define n^b = a \E^\. Then 


ICc 


(-1)1^1+^ IGd (g)(® ee ^A u(e) )Zl E »\ 

A 4 

( (^) ®e&E^K{e))Qc, (Bl) 

44 7^^ 


where Q c = i n * b X c if n^b is even and Q c = 
—i 1+n * b (—l) Sa c< x Y c if n^b is odd. This gives complete 
characterization of stabilizer generators, i.e., Q c = ±X C 
or Q c = ±Y C and the exact sign can be determined. This 
is essential in checking the incompatibility condition. 

The above discussion also justifies the assignment of 
the self-loop associated with the graph staet |Gq): 


n y mod 2 = 1, for T = x, 

n x mod 2 = 1, for T = y , 

n y mod 2 = 1, for T = z. 

Note that the stabilizer operators are not always in 
the canonical form in which K, c \ c = X c , i.e., they can be 
±X C or ±Y C , but those non-central operators are always 
Z. But it is easy to find rotations (around logical 0 -axis) 
to make them canonical. 


a c 

Z 

X 

y 

b 

X 

z 

Z 


V 

y 

X 


TABLE II. The choice of b and a^b- 


Moreover, as shown in the next Section, in terms of 
logical X, the measurement operator O c = (— l)\ Vc \X c , 
where the \V C \ denotes the number of sites in the domain 
c, so that the effective measurement by K on all sites of 
a domain gives rise to a projection P c = (l + O c )/2. Tak¬ 
ing into account of the probability it corresponds to an 
operator P c = (1 + O c )/ 2^1 (see below). If K\ c = ±X C , 
then the induced measurement by K a corresponds to X 
measurement in the canonical basis of the standard graph 
state formulation. If X c \ c = ±V C , then the induced mea¬ 
surement by K a corresponds to Y measurement in the 
canonical basis. Thus whether the induced measurement 
by K a on a domain is in the canonical basis either X or 
Y can be easily determined by whether n^b is even or 
odd. 


Appendix C: Proof of the weight formula 

Let us mention first the following fact that (b is chosen 
according to a c in Table 0. 

(G 0 1Orest = 0, (Cl) 

if I c is a strict subset of virtual qubits in any domain C 
(i.e. \I C \ < 4|C|) and ab is chosen according to Table [TT] 
(O r est denotes operators not in the support of domain 
C ). This can easily be proved by the fact that one can 
choose a stabilizer S :jq = X 3 X q a\ : !}al' 1 } (see Table]!]), where 
j £ I c and q £ C but q ^ / c , so that (<8>ie/ c (°|^) and Sj q 
anticommutes. Hence, 

(Go|Ore S t(® ie / c ai ll )|Go) 

= (Go|Orest(®ie/c CT I^)‘S'ifc|Oo) 

= — (Go|5 , jfcOrest(®ie/c Cr I^)|Oo) 

= -(Go|Orest(®<6/ c 4 1 )|G 0 ), 

showing that the expectation value is identically zero. 

Let us also note the following useful relation regarding 
to the 4-qubit GHZ associated with the corresponding 
POVM outcome K a , 


|GHZ")(GHZ"| =n ( 


(1- 


[i?;l] [i>;2] [v;3] [v;4] 


n n 


(C2) 

where n a (a = x,y,z) is a projection to a two- 
dimensional subspace, equivalently an identity opera¬ 
tor on the code subspace and can be safely omitted 
when acting on the graph state |G 0 ). Specificially, 
IBr = |+ + + +}(+ + + + | + | )( |, 
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YVy = + | - I, 

and n z = |0000)(0000| + 11 111)<11111. The label b a de¬ 
notes the corresponding type b if a c = a; see Table [TT| 
For a given domain (with a given type a), the POVM 
outcome on any site in the domain can be either F a or 
K a . Regarding the number nx of K outcomes, there are 
two scenarios: (i) uk is less than the total number \V C \ 
of sites in that domain C; (ii) uk = \Vc\- 

For case (i), the effect of all those K in terms of the 
probability distribution (or the weight formula) is to mul¬ 
tiply a factor of 2~ nK , i.e., (using J £ C to denotes the 
set of those sites with K ) 


is 


P({F,K}) 


0 f {v Af { v ( \ 0 k { A\k {iu i\\ 

sr T ol(v) ol(v) Zr T ac(w) ot(w)' 
v£Jf w€lJk 




I Jk\ 


vGJf 


]?( v ) 

ot(v ) a:(?;) 


0 

w£j K 


p{w)^ p{w)] p{w) p(w ) 

ol(w) ^ct(w) v a(w) a(w) 




I Jk\ 


0 

veC 


F y\ 

a{v) 


w£Jk 


|ghz; w )(ghz; w | 


F 


( U ) 


U^C 


a(u) 


W- 


(Go|O res t (0* eJ |GHZ Q(fl) ><GHZ a(ll) |) |G 0 


(! - A 


— (G 0 1Orest 0-ueJ 

= 2-^(G 0 |O rest |G 0 ), 


|G 0 ) 


In the second equality we have used the fact that K a = 
•\/3/2 K a F a , and in the third equality we have combined 
all F' s and written explicitly K a 's in terms of the four- 
qubit GHZ projectors. 

Now we know from Ref. 20 that 


where O res t denotes operators not in the support of do¬ 
main C, and we have used 



(C3) 


(G 0 |O 


rest v9v€lJ 


( >; 1 ]>;2] fT h ; 3] £T M] )| G o) = 0 . 


For case (ii), when we expand all the 2' v °' 
terms in 0„ 6 c(l — the only 

two nonvanishing contributions are 1/2^1 and 
(_l)l v d(04IG cr H)/2 | Vc | = (_l)IKIx c /2l Ve l. In terms 
of logical X, the effect of all K is equivalent to 
P c = (1 + O c )/ 2 l v tl ; where the O c = (— l)\ Vc \X c . That is 

(Go|O res t (® we c|GHZ a( . i; ^)(GHZ a(i) ^ |G 0 ) 


(1- 


— (Go I Orest 0i?eC 


= 2-1^1 (Go|Orest(l + O c )|Go). 


[v,l] [v-2] [v,3\ [v-A] 

U U Oi Ur U 1 


|G 0 ) 


where |Go) is an encoded graph state whose graph Go 
is specified by the POVM elements {F}, V is the set 
of domains of same-outcome POVM measurements, and 
£ is the set of inter-domain e dges (before the modulo-2 
operation) [21]. The formula (C3) was originally stated 


for the honeycomb lattice, but holds for all bipartite lat¬ 
tices. (For non-bipartite lattices, an additional condi¬ 
tion needs to be imposed relating to geometric frustra¬ 
tion [31j . Namely, if any domain contains a cycle with 
odd number of sites, such {F} will not appear.) Com¬ 
bining the above two results we find that 


p({ p , K }) = M 


\£\-\V\+\Jk\ 


x(G 0 | 


v&Jk 


|ghz- (w) ><ghz- 


( V ) 


I Go). (C4) 


Here we also see that the effect of all K in domain C 
is to measurement the logical qubit C in the logical X , 
followed by a post-selection of the result corresponding 
to either positive (if |V C | is even) or negative (if |V C | is 
oddd) eigenvalue of A'. 

With the above preparation, we can move on to the 
proof. Now consider a spin-2 AKLT state on a bi- 
colorable lattice C (generalization to non-bicolorable lat¬ 
tices is possible), and POVM elements F a and I\ a (a = 
x,y,z). Denote by Jf C £ the set of sites where the 
POVM outcome is of F-type and by Jk = £\Jf the set 
of sites where POVM outcome is of F-type. We should, 
strictly speaking, use a{v) to denote the type of x, y, z at 
site v. When there is no confusion, we simply write a. 


Using Eq. (C2) and the results in the beginning of the 
section, we know that for those GHZ-projections in a do¬ 
main such that their number is less than the total num¬ 
ber of sites in the domain, i.e., case (i) discussed above, 
their contribution is to mulitiply by a factor 2~ nK . For 
those such that the two numbers are equal, i.e., case (ii), 
these GHZ-projections (in a domain) can be replaced by 
P c = (1 + O c )/ 2l v °l, where the O c = (—l)l Vc lA c , and c 
labels the domain. Thus, 


\ \E\-\V\+2\J K \ 

p({F,K}) = \co \ 2 (^J 

x(G 0 1 0 (Ic + O c )\Go), (C5) 


Proof of Lemma 1. For simplicity let us denote the 
AKLT state by | ijj) below. The probability p({F, K}) for 
obtaining POVM measurements {F, K} described above 


where we use Dk to label the domains that contain the 
same number of I\ operators as the total number of in¬ 
ternal sites. 
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Next we demonstrate the first part of the Lemma. As¬ 
sume that, for some subset Q G Dk , the observable 
— 0 c g Q O c G iS(|Go))- Then, 

(G 0 | 0 {I/i + O^Go) 

= (Go| 0 (/„ + O v ) 0 {I)i + 0)i) (- 0 O c ) |G 0 ) 

V(zlD k\Q \ cGQ / 

= — (Go| 0 {Iy + 0^) 0 (O c + 7 C )|Go) 

v£D k \Q c(zQ 

= — (Go| 0 (d/i + O m )|Go) = 0. 

In the third line we have used the fact Oj) = J M . Let 
us also note that being product of Pauli operators, 
either commutes or anticommutes with another product 
of Pauli operators. 

Next, we demonstrate the second part of the Lemma, 
i.e., finding p({F, K}) when it is not identically zero. 
Consider a subset of domains Q C D K . If O m ^ 
±<S(|Go)), then (Go| 0^eQ O^IGo) = 0 (note that /r is an 
index for the domain, not an index for the site). Fur¬ 
thermore, if the incompatibility condition is not satis¬ 
fied, then 0 /i6 q 0)1 G ±<S(|G 0 )) implies that 0 (Ue Q O M G 
5(|Go)), and therefore (Go| 0 m gq O^IGo) = 1. We now 
exapnd the projector 0 ce d k {I c + O c ) in the matrix ele¬ 
ment, 

(Go| 0 {I c + O c )\G 0 ) (C6) 

cGD k 

= (G 0 | V 0 Ofi\G 0 ) = |M|, (C7) 

QCD K » eQ 

where the set M is defined as M = {0{Q) = 
®^qO w \Q C Dk and 0{Q) G <S(|G 0 ))}. Actually M 
has the following equivalent formulation which will turn 
out to be useful, 

M = {0{Q)= 0 OJQ C D k 
Q 

and [0{Q) 1 5] = 0, V5 1 G <S(|G 0 ))}. (C8) 


Using this latter characterization of M, we now turn to 
the counting for \M\. We describe every subset Q of Dk 
by its characteristic vector q, defined as follows: if p G Q 
then Q)i = 1, or if /i ^ Q, then = 0. Furthermore we 
define a binary-valued matrix H of dimension |U| x \Dk\ 
(where \V\ denotes total number of domains), whose en¬ 
tries are 


H, JV = 0, if [AC Mi O„] =0, 
H) 1V = 1, if {^,0^ = 0, 


where fi G V (the set of all domains) and v G Dk 
(the set of those domains with equal number of K's and 
sites). Then for any Q C Dk, 0{Q) G M if and only if 
-f/qmod2 = 0. Therefore, 

|M| =2 dim ( ker(ff) ). (C9) 

Putting everything into the expression for p({F, K}) we 
obtain the equation ([9]), 

/i \ |£|-|V|+2|Jir|-dim(ker(ff)) 

p{{F,K}) = \c 0 \ 2 3 4 5 6 7 

and the Lemma is proved. 

We remark that checking the kernel of a binary ma¬ 
trix can be done via, e.g., the Gauss elimination method; 
see e.g. j60j . Furthermore, to check the incompatibility 
condition it is sufficient to check the products of as¬ 
sociated with all basis vectors q’s in the kernel. If none 
of them satisifies it, then the incompatibility condition is 
not satisified. 
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